Coherent Development of Neutrino Flavor in the Supernova Environment 
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We calculate coherent neutrino and antineutrino flavor transformation in the supernova envi- 
ronment, for the first time including self- consistent coupling of intersecting neutrino/antineutrino 
trajectories. For neutrino mass-squared difference |<5m 2 | = 3 x 10 _3 eV 2 we find that in the nor- 
mal (inverted) mass hierarchy the more tangentially-propagating (radially-propagating) neutrinos 
and antineutrinos can initiate collective, simultaneous medium-enhanced flavor conversion of these 
particles across broad ranges of energy and propagation direction. Accompanying alterations in 
neutrino/antineutrino energy spectra and fluxes could affect supernova nucleosynthesis and the ex- 
pected neutrino signal. 

PACS numbers: 14.60.Pq, 97.60.Bw 
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In this letter we present the first self-consistent solu- 
tion to a long standing problem in following coherent fla- 
vor inter-conversion among neutrinos and antineutrinos 
in the region above the hot proto- neutron star subsequent 
to the supernova explosion [1]. The problem is that fla- 
vor histories on intersecting neutrino/antineutrino world 
lines can be coupled by neutral-current forward exchange 
scattering [2] . Many studies neglecting this aspect of fla- 
vor transformation 0, 0, 0, 0, 0, BS 03 nevertheless 
have shown that flavor conversion in the neutrino and 
antineutrino fields above the proto-neutron star could 
be important in understanding the supernova explosion 
mechanism J3|, |j. |5l and the origin of heavy r-process 

nuclei 0,0, 0,0, mm- 

Inelastic scattering processes and associated de- 
coherence may dominate neutrino flavor development in 
the proto-neutron star core and in the region near the 
neutrino sphere, necessitating a full quantum kinetic ap- 
proach there EH- By contrast, in the hot bubble 
(a high-entropy region that develops above the neutrino 
sphere at time post-core-bounce tpB > 3 s), where the 
r-process elements may be made, neutrinos and antineu- 




FIG. 1: Flavor evolution of neutrino v on trajectory desig- 
nated by angle $o relative to the neutron star surface normal 
is coupled with the flavor development of all neutrinos v' on 
intersecting trajectories. 



trinos for the most part propagate coherently. In this 
limit we can model the evolution of flavor along a single 
neutrino's trajectory with a mean field jl3l ]. Schrodinger- 
like equation, which in the effective 2x2 mixing channel 
is 



dt 



Hi 



(1) 



The effective neutrino flavor evolution Hamiltonian can 
be expressed in the flavor basis as 



-A cos 2(9 + A + B A sin 26 + B er 
A sin 2(9 + B* A cos 2(9 - A - B 



(2) 



The flavor evolution of an antineutrino is determined sim- 



ilarly but with A 
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-B and B e 
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In these expressions t is an Affine parameter along the 
neutrino's world line, A = Sm 2 /2E U , where E v is neu- 
trino energy, and 5m 2 = m 2 — m 2 is the difference of 
the squares of the relevant vacuum neutrino mass eigen- 
values. We focus on the atmospheric mass-squared dif- 
ference Sm 2 tm because it will give flavor transformation 
deeper in the supernova envelope than will the solar scale. 
Therefore, we set 5m 2 = ±3 x 10 _3 eV 2 , where the plus 
(minus) sign is for the normal (inverted) mass hierarchy. 
In Eq. 0, is the effective 2 x 2 vacuum mixing angle. 
Flavor transformation in the v e ^ z^ r and v e ^ 
channels is most important in supernovae because there 
may be disparities in energy spectra and fluxes among 
the neutrino flavors, and because v e and v e play a promi- 
nent role in setting composition and in the prospects for 
signal detection. In these channels ~ #13 for Sml tm . 
Experiment suggests sin 2 2#i3 < 0.1 and in our nu- 
merical calculations we take = 0.1. The flavor diago- 
nal potentials in Eq. ([2]) are A, from ^-electron forward 
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FIG. 2: (Color online) Plots of energy- averaged survival probabilities (Pw) for v e (left panels) and v e (right panels) as functions 
of radius r for the normal (upper panels) and inverted (lower panels) neutrino mass hierarchies, respectively. The solid and 
dashed lines give average survival probabilities along the trajectories with cos^o = 1 and cos^o = 0, respectively, as computed 
in the multi-angle simulations. The dotted lines give the average survival probabilities computed in the single-angle simulations. 



scattering (determined by the matter density profile in 
the hot bubble and B, from neutrino- neutrino for- 
ward scattering. The flavor off-diagonal potential B er 
similarly stems from neutrino-neutrino forward exchange 
scattering [l5| • Eq. (pQ) is nonlinear in that the potentials 
B and B er depend on the amplitudes a ea and a ra for a 
neutrino with initial flavor state a = e, r to be either elec- 
tron or tau flavor, respectively. However, the true com- 
plexity of this problem arises from quantum mechanical 
and geometrical coupling of neutrino / antineutrino flavor 
histories as illustrated in Fig. [TJ B and B er help de- 
termine the flavor development at point Q on neutrino 
z/s world line, but these potentials depend on a coherent 
sum over all neutrinos v' passing through Q. Here v T 
designates the relevant combination of mu and tau flavor 
neutrinos assuming these species are maximally mixed in 
vacuum and in the supernova medium [lfj. In our ex- 
ample numerical calculations we have taken the initial 
neutrino and antineutrino energy spectra to be of Fermi- 
Dirac form, with degeneracy parameter r\ v = 3 and av- 
erage energies (E Ue ) = HMeV, (E„ e ) = 16MeV, and 
(E Ut ^ t ) = 25MeV and we take the energy luminosity 



for each neutrino species to be L v = 10 51 ergs 1 . 

Though neutrinos and antineutrinos are emitted from 
the neutrino sphere (radius r = 11 km) in a thermal, in- 
coherent manner, our simulations show that large scale 
coherent and collective flavor transformation develops 
with increasing radius r. This behavior is driven in 
part by progressive forward scattering-induced coupling 
of flavor evolution on intersecting neutrino and antineu- 
trino world lines. Very far from the neutron star sur- 
face, where neutrino densities are low, B and B er are 
negligible .and ordinary Mikheyev-Smirnov-Wolfenstein 
(MSW) [I?], EH, EH flavor evolution or vacuum oscilla- 
tions obtain. However, at smaller r, B and B er can 
dominate and flavor evolution can be very different from 
MSW. For the normal mass hierarchy, our simulations 
show that the more tangentially-propagating neutrinos 
and antineutrinos, which experience the largest B and 
B er potentials because of the intersection angle depen- 
dence in the weak current @, 0] , are the first to experi- 
ence significant flavor transformation for a broad range 
of energies. This sets in, e.g., at r « 80 km if the entropy- 
per-baryon in the hot bubble is S = 140 in units of Boltz- 
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FIG. 3: (Color online) Plots of survival probabilities P uu for neutrinos (left panels) and antineutrinos (right panels) as functions 
of both neutrino energy sij-u and emission angle i9o at radius r = 225 km. The upper panels employ a normal neutrino mass 
hierarchy, and the lower panels employ an inverted neutrino mass hierarchy. 



mann's constant k^. This simultaneous conversion of u a 
and u a quickly spreads to all neutrino/antineutrino tra- 
jectories, leading to coherent, collective flavor oscillations 
of the entire neutrino/antineutrino field. For the inverted 
mass hierarchy, the opposite is true: radially-propagating 
neutrinos and antineutrinos transform first. 

These features of flavor development can be seen in 
Fig. [2j The survival probability at location t along a 
given neutrino's world line is, e.g., for a neutrino which 
is initially electron flavor, P UeUe (t,$o,E u ) = \a ee (t) | 2 . 
In Fig. [2] we show the energy-spectrum- averaged survival 
probabilities (Puu) for v e and v e as functions of r for 
both the normal and the inverted neutrino mass hierar- 
chy cases. Here the energy averages are over the initial 
energy spectra for each flavor. It is clear that flavor evo- 
lution along different trajectories can be different, yet it 
is also evident that neutrinos and antineutrinos can un- 
dergo simultaneous, significant medium-enhanced flavor 
conversion. Our simulations show that this conversion 
can take place over broad ranges of neutrino and an- 



tineutrino energy. We have also performed simulations 
using the single-angle approximation widely adopted in 
the literature. These give results qualitatively similar to 
our multi-angle calculations, as shown in Fig. [2l The 
collective neutrino flavor transformation observed in our 
simulations is not the "synchronized" mode described in 
Ref. [7]. In the normal mass hierarchy case, neutrinos 
or antineutrinos in the synchronized mode undergo one- 
time transformation in the same way as does a neutrino 
with energy p syYl c • There would be little synchronized 
flavor transformation in the inverted neutrino mass hier- 
archy. 

The collective neutrino flavor transformation evident 
in Figjp] is likely of the "bi-polar" type as described in 
Ref. [20]. In this mode, neutrinos and antineutrinos ex- 
perience in-phase, collective, semi-periodic flavor oscilla- 
tions, even for the inverted mass hierarchy. This behavior 
was first observed in numerical simulations of neutrino 
flavor transformation in the early universe 21, [22[. It 
has been argued [20] that neutrinos and antineutrinos 
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could undergo collective flavor transformation in the bi- 
polar mode for the inverted mass hierarchy case even in 
the presence of a dominant matter potential. This con- 
jecture seems to be supported by our simulations. 

Using the single- angle approximation, we have stud- 
ied the relation between rx, the radius where large-scale 
flavor transformation first occurs, and the entropy per 
baryon S and neutrino luminosity L v . With larger S 
the matter density profile is more condensed toward the 
neutrino sphere. We find that, for the normal neutrino 
mass hierarchy case, rx decreases substantially if S is 
increased from 140 to 250. However, the value of rx de- 
creases only slightly with the same change in S for the 
inverted neutrino mass hierarchy case. The non-linear 
effect of the neutrino-neutrino forward scattering poten- 
tials is enhanced with higher L v . For the normal neutrino 
mass hierarchy case, rx decreases monotonically with in- 
creasing L V1 and approaches the radius where neutrinos 
and antineutrinos start collective flavor transformation in 
the synchronized mode. For the inverted neutrino mass 
hierarchy, rx increases monotonically with increasing L v . 

In Fig. [3] we show survival probabilities P vv as func- 
tions of both neutrino energy E v and emission angle $o 
at radius r = 225 km. Large-scale quantum interference 
stemming from coupling of flavor histories is evident. The 
horizontal "fringes" in Fig. [3] are generated during neu- 
trino/antineutrino background-dominated collective fla- 
vor transformation. In this regime E v is unimportant, 
but $o is crucial. This is because the neutrino-neutrino 
forward scattering potentials are energy-blind but have 
a strong angle dependence. The vertical "fringes" are 
generated at larger r where the neutrino background is 
weaker and, therefore, (non-adiabatic) MSW transforma- 
tion dominates flavor evolution. In this case $o is unim- 
portant and E v is crucial. This is because at large r the 
neutrino beams are almost parallel to each other, and the 
MSW effect is energy dependent. Our simulations show 
sharp, vertical transition regions at E v c± 10 MeV for 
neutrinos in both the normal and inverted neutrino mass 
hierarchies. This feature is likely a result of break-down 
of collective transformation [l[. 

Although the establishment and breakdown of collec- 
tivity in v and v flavor transformation remains an open 
issue for the general supernova environment, our particu- 
lar simulations are based on reasonable assumptions and 
our results are robust with these assumptions. These 
simulations suggest that, with Sm1 tm , large-scale collec- 
tive flavor transformation in the bi-polar mode can occur 
deep in the supernova envelope. In broad brush, across 
all $o mostly lower energy v e and few v e are transformed 
in the normal mass hierarchy. The opposite is true in 
the inverted mass hierarchy. However, as is evident in 
Fig. [3] at r = 225 km, in either mass hierarchy, survival 
probability can show significant $o dependence. The v e 
and v e energy spectra and angular distributions could 
be quite different from those in zero-neutrino-mass su- 



pernova models. This could affect the expected neutrino 
signal and conceivably affect neutron-proton interconver- 
sion rates and the prospects for r-process nucleosynthesis 
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